function tech = technology(t,alpha, hc, I, PS, H, TFP_param)

TFP = exp(  TFP_param(1) + TFP_param(2)*(t) + alpha(6).*PS );
%If P=1;
theta_tp1_PS1 = ( TFP.*[ alpha(2).*hc.^(alpha(4)) +  (1-alpha(2)).*[ alpha(3).*H.^(alpha(1)) + (1-alpha(3)).*I.^(alpha(1)) ].^((alpha(4))/alpha(1)) ].^(alpha(5)/alpha(4)) ).*(PS==1) ;
%If P=0;
theta_tp1_PS0 = ( TFP.*[ alpha(8).*hc.^(alpha(10)) +  (1-alpha(8)).*[ alpha(9).*H.^(alpha(7)) + (1-alpha(9)).*I.^(alpha(7)) ].^((alpha(10))/alpha(7)) ].^(alpha(11)/alpha(10)) ).*(PS==0) ;

theta_tp1 =  theta_tp1_PS0 + theta_tp1_PS1;

tech = theta_tp1;